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ABSTRACT 

'(Baroclinic eddies with diameters of about 10 km and 
current speeds of about 50 cm/s have been widely observed in the 
central Arctic Ocean north of Alaska and Canada. The possible 
origin of these eddies through an instability of the mean baro- 
clinic flow is investigated using an ocean model with exponential 
profiles of mean shear and VSis&ia frequency. The model includes 
Ekman pumping at a rigid bottom and at either a free or rigid 
upper surface. The central Arctic Ocean where the eddies were 
found is baroclinically stable with no possibility of eddy pro¬ 
duction. If the eddies are spawned by this mechanism, they must 
be formed at a site far from where they were observed. On the 
periphery of the Arctic Ocean north of Alaska the combination of 
greater current shear, shallower depth and lack of ice cover 
leads to unstable conditions and the eddies apparently originate 
in that region. The instability theory predicts maximum velocity 
at the surface instead of below the surface as observed. Appar¬ 
ently after formation in open water the eddies are advected 
beneath the ice cover and dissipate the momentum of their upper 
layers against the ice. This is demonstrated by calculations for 
the diffusion of vorticity against the ice in the case of an 
initial exponential profile. A subsurface maximum then develops 
which resembles the observed profiles. 
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1. Introduction. 


The temperature, salinity and velocity fields of the oceans 
have a mesoscale variability on scales of tens to hundreds of kilo¬ 
meters horizontally and on the order of months in time. These 
mesoscale features frequently take the form of approximately circu¬ 
lar baroclinic eddies with a nearly geostrophic force balance. It 
is likely that mesoscale currents play an important part in driving 
the large-scale ocean circulation and in transporting heat as well 
as chemical and biological parameters on a global scale. A review 
of the literature of mesoscale currents (McWilliams, 1979) shows 
that these features have been widely observed and appear to be 
nearly ubiquitous in the world ocean. 

Mesoscale eddies in the Arctic Ocean were reported by Hunkins 
(1974) and Newton et al (1974) based on research carried out from 
drifting ice stations. During the AIDJEX (Arctic Ice Dynamics 
Joint Experiment) project in 1975 and 1976, a much more extensive 
survey was made and large numbers of arctic mesoscale eddies were 
identified (Hunkins, 1980; Manley, 1981). Profiles of temperature, 
salinity and current were made on a daily basis from four manned 
AIDJEX ice camps spaced 100 km apart drifting in the central Arctic 
Ocean. The arctic eddies are strongly baroclinic with a subsurface 
velocity maximum in the steep density gradient between 30 and 200 m. 
Vertical profiles of oceanographic parameters through an arctic 
eddy and in adjacent waters outside the eddy are shown in Figs. 1 
and 2. The maximum current speed of 58 cm/s in this particular 
eddy contrasts with the geostrophic mean currents in this region 
which are only 5 to 10 cm/s. During their 14-month drift which 
covered a total track length of 8,140 km, the four AIDJEX camps en¬ 
countered 143 eddies. The location of these eddies is shown in Fig. 3. 








One of the eddies was recrossed after 5 h months, demonstrating that 
they persist for at least that length of time. There is a peak in 
mean vertical kinetic energy profiles at 120 m coinciding with the 
velocity maximum. The kinetic energy of the central Arctic Ocean 
resides almost entirely in the eddies with the mean currents ac¬ 
counting for less than one per cent. Although detailed surveys of 
individual eddies are lacking, most of them seem to have a diameter 
of about 10 km. On that basis, it is estimated that eddies cover 
between 10 and 20 per cent of the area of the western Arctic Ocean. 
Their sense of rotation is anticyclonic or clockwise. Water prop¬ 
erties within the arctic eddies differ from those without and 
suggest that these features originate in the shallow waters of the 
Chukchi Sea. 

A number of possible generating mechanisms for these eddies 
have been suggested including the stress of wind and drifting ice, 
brine convection accompanying sea ice growth, and an instability 
of the mean baroclinic current. Of these, the baroclinic insta¬ 
bility hypothesis seems most plausible and is the only one examined 
here in detail. 

It is generally accepted that the migrating cyclonic systems 
of the atmosphere as well as the turbulent flow in certain labora¬ 
tory experiments with rotating fluids result from the instability 
to small disturbances of the mean stratified shear flow. This 
identification of baroclinic instability as a mechanism for gener¬ 
ating eddies in the atmosphere and laboratory has produced a body 
of theory which has been summarized and reviewed by Charney (1973) 
and Kuo (1973). It has also been suggested that mesoscale eddies 
observed in the ocean may be initiated by the same instability 
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mechanism (Orlanski and Cox, 1973; Gill et al, 1974). 

One of the simplest theoretical demonstrations of baroclinic 
instability is that by Eady (1949) who showed that for a basic 
state in which the vertical shear, vaisSla frequency and Coriolis 
parameter are all constant, certain intermediate wavelengths are 
unstable and tend to grow. Since vertical shear and VaisalS fre¬ 
quency are not actually constant in the ocean, models with more 
realistic structure are desirable for comparison with data. 

Eady's theory was generalized by Williams (1974) to include mean 
vertical profiles which vary according to a power law. He further 
included frictional effects by introducing Ekman pumping at rigid 
upper and lower boundaries (Williams and Robinson, 1974). The 
results were successfully applied to interpret rotating tank ex¬ 
periments. For both Eady's and Williams' problems, solutions are 
in terms of hyperbolic functions. 

In the present paper Eady's problem is extended to include 
mean shear and vaisaia frequency which vary exponentially with 
depth. Analytical solutions to the problem are in the form of mod¬ 
ified Bessel functions. Friction is included in the form of Ekman 
pumping at either rigid or free surfaces. A model similar to that 
developed, in this paper has been described by Hart and Killworth 
(1976) who also applied their results to the Arctic Ocean. The 
model used here differs by including free-surface Ekman pumping and 
a discussion of vertical wave structure. 

Observed eddies in the Atlantic Ocean are in the neighborhood 
of 200 km in diameter and apparently have their amplitude maximum at 
the surface. The arctic eddies are only about 10 km in diameter and 
have their maximum amplitude in the pycnocline well below the 
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surface. These differences between eddies in the Atlantic and 
Arctic Oceans suggested an investigation of factors influencing the 
horizontal scale and vertical profile of the disturbance. The 
shape of the mean shear and stability profiles as well as the pres¬ 
ence of the ice cover on the Arctic Ocean are expected to be impor¬ 
tant factors. 


2. The Exponential Eady Problem . 

The analysis is based on the quasi-geostrophic potential vor- 
ticity equation in a form used extensively in oceanography and 
meteorology. For an ocean with a mean state varying only with 
depth and which is disturbed by small perturbations in pressure and 
density we have 

a " d = f. (Z.) + /**' (^ fc ) 


where the perturbations, p* and , are much smaller than the mean 
values. If there is a geostrophic balance so that 


U = 


AT = 


_ L IV 
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where is tii e geostrophic stream function. 
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then the quasi-geostrophic equation may be written 


T> 

Dt 


v"r +■ 





o 

( 2 . 1 ) 


where directions are x, east; y, north; z, vertically upward; and 
the north and east velocity components are u and v. The Coriolis 
parameter, f, has its usual meaning. Other symbols are defined as 
follows: 


Total derivative, 


5- 



LA \ — 4- N~ 



Horizontal Laplacian, 


vaisala frequency, 


Beta parameter. 




a/f 

i ~z- 


The quantity in brackets is the potential vorticity, composed 
of the relative vorticity and the stretching affect of density 
gradients on vortex lines. The final term of the equation repre¬ 
sents the effect of meridional motions on potential vorticity. For 
the derivation of this equation and discussion, see Pedlosky (1979). 

Assuming a stream function composed of steady and perturbation 


parts 



l/(z) 
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where the steady eastward flow varies only with depth, the linear¬ 
ized form of (2.1) becomes 


at: 




A£ ±.f£lw) 

ax J 



( 2 . 2 ) 

This equation may be separated by introducing a perturbation 
stream function in the form of a zonally-propagating wave, 

ikU-c-*-) 

f = "R,. F(z) si-w & 


to give the differential equation for vertical structure 




+ 


A i ft' AH) 
p - Tz. hr- 72 


(2.3) 


2 

If both N and dU/dz have the same functional form and 6 = 0, 

the final term in brackets vanishes and the singularity at U = c is 

circumvented. This has been called the generalized Eady problem by 

? 

Williams (1974) who used a power-law dependence for N" and dU/dz. 
Here we introduce an exponential depth dependence for those quanti¬ 


ties. 

Lot the basic density variation take an exponential form 

/?(z) - /f -P A/* (/ - 


where o s is the density value at the surface and Ap is the differ 
ence in density between the surface and great depths. 











The vaisaia frequency in this case becomes, 



c| o< A P 6^ _ ^ dd.A/° 4Z 

- ~*7 S 


to a good degree of approximation since, for the ocean, 






A / 5 = 1 Xj'*' 3 


and the scale depth of the pycnocline, 

o < = 10~ 3 m ' 1 


The value of the VaisSlS frequency at the surface is 

d < A 

t zlJ -— 

so that 




N 


x 



o< 2 

C 


The basic current takes the same exponential form, 

17 = 17 , 

These profiles have a resemblance to those in many regions of the 
ocean including the Arctic. It can be shown that the association of 
exponential stratification and exponential shear implies isopycnal 
surfaces which slope uniformly. Note that for these mean conditions 
the Richardson number increases exponentially with depth, 
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Introducing these exponential mean conditions into (2.2) produces 


the non-dimensional differential equation, 


cTF 


A x r— 

K F 

A z 1 

Jz. 


Circumflexes indicate 

parameters which have been 

of the pycnocline, z = 

= az, and by an internal Ro 

formation. 
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Introduction of a change of variable, 


F(z.) 



where 

A i/2^ 


t = 

2 X ^ 


transforms (2.4) into 

a standard form of 

the Besi 

<TP + 

_L iZ. - ( 

' + T 

A f 1 

t 1 
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(2.4) 


=< f 


m 
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(2.5) 


which has a solution in terms of modified Bessel functions of the 


= f\z, Cf) + ^ 


first order 





where A and B are arbitrary constants. Returning to the original 
dependent variable, the solution becomes 


FC-fO = Afl'Cf) + 

( 2 . 6 ) 


Boundary conditions at the surface and bottom are imposed by 
Ekman pumping as introduced by Charney and Eliassen (1949). Vor- 
ticity of the interior frictionless flow leads to horizontal mass 
divergence in the Ekman boundary layer and consequent vertical flow 
at the interface between the interior and the boundary layer. Thus 
frictional effects in the boundary layers are communicated to the 
interior as vertical velocities. The Ekman pumping boundary condi¬ 
tions imposed on the vertical velocity. 




are either 


_ JL 

H" Dt\^Z ) 
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for a rigid upper surface, 
such as an ice cover 

(2.7) 


or 


CM = - 


for a stress-free upper 
surface as in the case of 
an open ocean 


( 2 . 8 ) 







For a lower rigid surface the Ekman pumping condition is 


u = ^ Is 


The geostrophic vorticity is defined as 


r - v'f = - (k''+■»<]) f 


The scale depth of the frictional boundary layer is 

Dt = 

where ^ is a kinematic eddy viscosity. 
Substitution of (2.6) into the boundary conditions gives two 
equations in terms of the two unknown amplitudes, A and B. For rigid 
upper and lower boundaries, the equations are 

A [ C' -1) * /.(f.) - t. (f •) -IY it, (f .)] 

. $ [(£-/)* K.(f.) - K,(f.) - * ^t-)] = ° 


(2.9) 
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( 2 . 10 ) 


where 


r ZV 


a - 8/a- 

= 2.Vt 








while for a free upper surface (2.9) is replaced by 

A [(/-£) v Up) - Up) -Lwrx. (f.)] 

+ b[ (c-o v /^.(f.)-K.(f.) + 

where H is the dimensionless ocean depth, H , and the dimension¬ 
less phase velocity is defined by c = c/U Q . The Ekman pumping 
parameters are 

Y J>e Ns AXJ, 

for a rigid surface 

and 

W =K“<Y s/ ^i 7 , 

for a free surface. 

For other than trivial solutions, the determinant of the co¬ 
efficients, A and B, must vanish. This provides a dispersion equa¬ 
tion relating c and For real values of 9, c is generally complex 
with the real part, c^,, representing the phase speed and the imagi¬ 
nary part, c^, representing the growth or decay of the wave. 

Positive values of c^ indicate growing waves or instability. In 
this case the dispersion equation is a quadratic with two roots and 
the solutions consist of either a pair of neutral waves, both with 
c^ = 0, or one growing and one decaying wave. The roots, calculated 
with the aid of a digital computer, are plotted in Fig. 4 for H » 5 
and in Fig. 5 for H = 10. Three types of boundary conditions were 
used: frictionless, free-rigid and rigid-rigid. The phase speeds, 

c r , are plotted for both the growing and decaying waves. Note that 
for free-rigid boundary conditions the lower left-hand branch of 





the phase speed curve is associated with growing waves while for 
the rigid-rigid case it is the upper branch which corresponds to 
growing waves. The growth rate, ?ci, is plotted only for the 
growing waves. The examples all show a maximum growth rate at 
some intermediate wave number and a critical wave number beyond 
which no growth takes place. These results are qualitatively sim¬ 
ilar to results by others for Eady-type models. When the wave is 
unchanging in the y-direction, m = 0 and y = k to give the highest 
values of the wave number for the growth maximum and for the 
critical point beyond which growth ceases. The scaled water depth 
in Fig. 4, H = 5, corresponds generally with mid-ocean conditions 
where the scale depth of the pycnocline is about 1 km and the total 
water depth about 5 km. 

The highest growth rate is found for the case with no Ekman 
pumping where yc i * 0.134 at a wave number of 0.8 when H = 5. 

Beyond the critical wave number, ? = 1.18, no growth occurs and 
only neutral waves exist. When Ekman pumping is introduced the 
growth rate is reduced and the wave numbers of the maximum and 
critical points are shifted to lower values. For values of Y * 0.5 
and W =0.05, chosen to correspond roughly to oceanic conditions, the 
dimensionless growth rate is reduced to 0.068 for a free upper sur¬ 
face and to 0.02 for a rigid upper surface, with rigid bottom sur¬ 
faces in both cases. For a larger ratio of pycnocline depth to 
total depth, H = 10, growth rates as well as phase speeds are less 
(Fig. 5). 

The phase speeds of the decaying and growing waves coincide in 
the frictionless case, while the phase speeds of the neutral waves 
split into two separate branches. The presence of Ekman layers 
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splits the phase speeds for the growing and decaying waves also into 
two separate branches. 

Vertical profiles for the fastest growing wave are shown in 
Figs. 6, 7 and 8 for three different boundary conditions. The rela¬ 
tive perturbation stream function may be represented as 

t Ooi k ( x ' ct ) + i. sin k(x-tt) 

where the depth function, F, has been decomposed into real and 
imaginary parts. The relative magnitude and phase of the perturba¬ 
tion velocity are given by 

/ i. / x \ '/a- 

a\ = (a +b' ) 

and 

e * tan'' (B'/A') 

The perturbation velocity at the upper surface is greater than 
that on the lower surface in all of these cases, although in Eady's 
original problem the velocity profile was symmetrical about the 
middle level. The solutions are monotonic functions of depth in 
both problems and there are no maxima at intermediate depths either 
with or without Ekman layers. 

The negative phase which is plotted in the figures presents a 
lateral view of wave tilt in an east-west direction. The backward 
tilt of growing waves as the surface is approached is characteristic 
of unstable baroclinic waves. 

Linear perturbation theory predicts conditions only at the on¬ 
set of instability and does not describe later developments which 
occur as the unstable waves continue to amplify. However, it is 
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expected that the wave which initially grows fastest will continue 
to dominate the spectrum of unstable waves and that it will maintain 
the general character of its original profile. 

The theory presented here does not include the beta-effect as 
an approximation to the latitudinal variation of the Coriolis 
parameter. This neglect seems justifiable in the present study 
since it has been shown by Green (1960) in numerical investigations 
that although the inclusion of a finite beta-parameter in the Eady 
problem produces an additional weaker instability maximum at low 
wave numbers and also destabilizes the high wave numbers, it does 
not significantly alter either the wave number or the growth rate 
of the fastest-growing wave which remains essentially the same as 
in Eady's problem. 

3. Application to the Arctic Ocean . 

First we examine the possibility that the arctic eddies are 
generated in offshore locations near where they are observed. 

Examples of profiles for geostrophic current and vaisMla frequency 
plotted in Figs. 9 and 10 show that an exponential curve roughly 
fits conditions in the central Arctic Ocean. Appropriate constants 
based on these data are: 

A/. = /O ' 1 s -' v = a y 10 '*- brj-s ' 1 

° > s 

r) =/0vw - S X / o~ 3 m~‘ hi “ £ ooo k 

t J J J 

f - !.**■ *I°' H s ■' 


The Ekman pumping parameters are then 

Y- r and W r O.Z5- 
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The exponential model is stable for Y > 1.5 in the case of a 
rigid top representing the permanent ice cover for H = 10. Hence 
the central Arctic Ocean is stable to small disturbances and the 
eddies cannot originate there by this process. If we choose an 
improbably low value for the Fkman pumping parameter of Y = 0.5, the 
flow is unstable but the growth rate 

X Ci 

corresponding to an e-folding time of 612 days, is much too long for 
significant growth in this situation. The half wavelength for the 
fastest-growing wave in this case would be 

Va = 

which is an order of magnitude larger than the observed diameter. 

If the arctic eddies are formed in the Arctic Ocean by baro- 
clinic instability, they must originate at some other location 
where conditions are less stable. At the edge of the Alaskan Con¬ 
tinental shelf, higher mean shear, shallower depth and lack of ice 
cover during late summer all enhance instability. An oceanographic 
section made by the USS Staten Island in September 1959 across the 
Alaskan shelf north of Point Barrow (U. S. Naval Oceanographic 
Office, 1963) provided guidance in selecting parameters for this 
region: N 3 ■ 10 “ 2 S ~ 1 , U g = 10”^m.s“^, D Q * 10m, a = 10 -_ m - ^, 

H * 1,000m (Y = 0.5 and W = 0.05). For free upper and rigid lower 
boundries, the flow is unstable with a scaled maximum growth rate 
Of 0.055 at a scaled wave number of 0.6 corresponding to a dimen- 
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sional growth rate of 15 days with a half wavelenght of 37 km. 

This rapid growth rate indicates that mesoscale eddies should form 
in this area. The dimensions of the predicted growing disturbance 
are roughly comparable to observed eddy diameters. Eddies thus 
formed would be advected by large-scale currents to the sites 
where they were observed in the deep ocean away from shelf areas. 

On the other side of the Arctic Ocean in the area northwest of 
Svalbard there is another region of strong shear which can also be 
examined for instability. Warm water from the Atlantic flows north¬ 
ward into the Arctic Ocean on the eastern side of the strait be¬ 
tween Greenland and Svalbard while cold arctic surface water leaves 
on the western side. Values for the critical parameters in that 
area may be assigned on the basis of a survey by Newton (1977): 

N g = 10" 2 s -1 , U s = 0.065m«s - ^, D e « 7m, a = 1.25xl0 -2 m~^ and 
H = 800m so that y = 0.5 and w = 0.05. The flow is unstable with an 
e-folding time of 18 days for the fastest growing half wavelength 
of 30 km. These results are similar to those in the Alaskan shelf 
region but it is not likely that the observed eddies originate so 
far away. Oceanographic surveys are being planned in the eastern 
Arctic Ocean which should provide information on eddy population 
north of Greenland and Svalbard. Eddies there would be involved in 
the exchange of heat and salt between the Atlantic and Arctic Oceans 
and may be a factor in climatic change. 

4. Vorticity diffusion and the subsurface current . 

Current speeds in the disturbances generated by baroclinic in¬ 
stability have their maximum at the top pf the model which corres¬ 
ponds to the base of the Ekman layer at a depth of about 10 m. Some 
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explanation is required for the fact that the observed velocity 
maximum is much deeper, averaging 120 m. If the eddies are indeed 
generated along the edge of the Alaskan shelf during summer when 
open water prevails there, then the current maximum will be near 
the surface at the base of the Ekman layer. If these eddies are 
then carried by large-scale currents into the central ice-covered 
ocean, the high surface velocities will decrease rapidly through 
friction against the ice leaving a subsurface maximum. This 
hypothesis can be considered quantitatively as a problem of upward 
vorticity diffusion from an eddy into a surface layer maintained 
at zero vorticity. 

First we need to examine the relation between the Ekman layer 
and the mixed layer which extends downward to a depth of about 
50 m. The mixed layer, which is maintained throughout most of the 
year by brine convection resulting from ice formation, will be 
spun down rapidly through friction against the overlying ice cover. 
Secondary circulations will cause an efficient Ekman spindown of 
the entire layer on a time scale of 

% s //*(2 

where H m is the mixed layer depth (Holton, 1972). Choosing repre¬ 
sentative values for the Arctic Ocean of H m = 40 m, f = 1.4 X 10 -4 s - * 
and 7 X 10~^m^s“^, we find a rapid response time of 16 hr. 

Below the mixed layer the strong stratification will suppress 
secondary circulation so that momentum and vorticity will be lost 
upward only by vertical eddy diffusion. Then the time scale will be 

X s K/» 
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where 11^ is taken as the distance from the base of the mixed layer to 

the subsurface velocity peak, about 100 m, and ^ must have a re- 

-4 2 

duced value of, say, 10 m /s so that 

% = 3 yr* 

Thus the response time of the mixed layer at the surface is negli¬ 
gibly small in comparison with that of the stratified layers below 
it. We may then assume that the frictional effect of the ice is 
transmitted directly to the base of the mixed layer. 

We examine the case of an initial vorticity profile with a 
maximum at the top, that is at the base of the mixed layer, where 
the vorticity remains zero due to friction and trace the changes in 
the profile with time as a subsurface maximum develops which resem¬ 
bles the observations. Since the horizontal and vertical structures 
are separable and since vertical vorticity is a horizontal deriva¬ 
tive of velocity, the vertical profile derived for vorticity will 
also apply to velocity. 

In the absence of vertical stretching and horizontal divergence, 
the linearized equation for vertical vorticity is a diffusion equa¬ 
tion, 

(4.1) 

where vorticity 

^ a*t~ _ «=) 
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In terms of a Fourier integral, (4.1) has a solution 


f- 


- o 


and 




—> 0 


■Jh*) 

.(z-z'/Aj it , 

c Jz 

cJ 

— aO 

(4.2) 

y is f(z') • 

The boundary conditions are 


JZ. - o for* t- c> 

d s 

yr — oO ■^V' dtt 


Although the vertical profile derived from the previous baroclinic 
instability problem has the form of a modified Bessel function, we 
approximate that shape here with an exponential. Assume an expo¬ 
nential profile which has a mirror image above to maintain zero 
vorticity at z = 0. , 

f(z!) = * *** 2> ° 

_ _ -z 



(z-l')Wtt 



-»z'«_ 

c. cPZ 
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or 


Z/z i) --j, -L 


»z +»* -* / -E z - c 2 + z. * *■ *)jy*»* 

t Jl 


.,Z+At«’/' -£*'-(Z-Zi, 

+ t ft 


(4.3) 


In the first integral of (4.3), let 


£lr (z+Z»t*) _ _ 
2 O t; 'A 


and in the second integral, let 
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Introducing the error function notation. 


e-r-p 


= £ 
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Jo 


we have 


7 (*» ■ 
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-e 


/ v- erf(- 


Z + lD-t* 
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/- £ ' 


!> / 2»t<K-Z 

'l 20>t)^ 


(4.4) 


Equation (4.4) may be evaluated in terms of normalized variables: 
for vorticity, 7// ; depth, az; and time, (\ ~t* ^ (Fig. 11) . 
Vorticity diffuses upward and downward so that a subsurface maximum 
develops which both decreases with time and descends to deeper 
levels. These theoretical profiles closely resemble the observed 
profiles, for example see Fig. 1. 

From the previous results on baroclinic instability we may 
assume that 1/a = 200 m. The choice of an eddy coefficient is more 
difficult but experience in steep density gradients of other oceans 
suggests a value of )) * 10”*mVs for an order of magnitude cal¬ 


culation. In the Arctic Ocean the velocity maximum occurs at 100 
to ISO m. The mixed layer depth of 50 m must be subtracted to cor¬ 
respond to the theoretical problem. Then az = 0.25 to 0.5. 






Referring to Fig. 11 it is seen that this corresponds to a dimen¬ 
sional decay time ranging from 46 to 185 days, indicating that the 
observed eddies had spent from 1 to 6 months beneath the ice. The 


decay time is directly proportional to choice of the eddy coeffi¬ 
cient value which is admittedly a poorly known quantity and this 
eddy age can only be considered suggestive. 

5. Summary and discussion . 

The vigorous mesoscale features which are so abundant in the 
Arctic Ocean require a source of energy for their generation and 
maintenance. Although wind energy has been suggested as a source, 
the lack of correlation between winds and mesoscale motion makes 
this an unlikely relationship. Furthermore the arctic eddies are 
much smaller in horizontal scale than synoptic weather systems 
although a similar size might be expected if the eddies were 
driven by winds. Another possible energy source which has been 
suggested is thermodynamic in nature and associated with the freez¬ 
ing process. During winter, areas of open water appear in the ice 
pack through either breakup or differential motion of ice floes. 

Tee grows quickly in these open areas at the extremely low air 
temperatures which prevail. As the ice grows, salt is expelled and 
' heavy brine sinks through the mixed layer causing intense local 

convection. 3ut it is difficult to understand how this shallow 
convection can produce the observed redistribution of density at 
the much deeper levels where eddies are found. Also there is here 
again a mismatch in scales between the sizes of source and eddy. 
Open areas where rapid freezing occurs are usually tens to hundreds 

of meters across while the eddies are about 10 kilometers in 
i 
i 
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diameter. A more plausible source for eddy kinetic energy than 
either of these suggestions is the available potential energy re¬ 
siding in the inclined density surfaces accompanying the mean 
geostrophic circulation of the Arctic Ocean. The mean vertical 
shear has been shown to be unstable in the Barrow Ov-ent north of 
Alaska but not in the central Arctic Ocean where eddies were ob¬ 
served. The small-amplitude theory employed predicts an eddy size 
close to that observed and it predicts a growth rate fast enough 
to indicate rapid formation of these features. Although the theory 
cannot describe the evolution of an eddy, it is assumed that the 
unstable waves continue to develop until meanders form which 
eventually cut off to become isolated eddies. 

The baroclinic instability theory used here applies to a ver¬ 
tical shear only and is appropriate to mid-ocean circulation where 
the horizontal shear may be considered negligible. It incorporates 
an exponential mean shear and stratification as well as Fkman 
pumping in an attempt to more closely fit actual oceanic conditions. 
The conclusion that the central Arctic Ocean is stable and that 
eddies do not originate there by this process seems quite definite. 
They must be formed elsewhere, apparently in the Barrow Current 
and application of the theory to this current confirms that it 
could be the source. 

There are two assumptions inherent in the theory which are not 
exactly met in the Barrow Current however. The Current is a rela¬ 
tively narrow jet which has significant horizontal shear and it 
occurs along the shelf edge where the topography is not flat. But 
it is not believed that the conclusions reached would be seriously 
altered by inclusion of a horizontal mean shear effects and depth 







variations. Although no quantitative modeling demonstrates this, 
some qualitative arguments make the omission of these effects ap¬ 
pear relatively unimportant. 

A mean horizontal shear may lead to instability and it is 
quite possible that the 3arrow Current is barotropically as well 
as baroclinically unstable. The combined barotropic and baroclinic 
instability problem is an inseparable one and solutions to it have 
been primarily sought by numerical methods. Barotropic instability 
could only be expected to enhance eddy production. Its effect on 
growth rate and eddy size in this case will require further model¬ 
ing in the future. 

3ottom topography such as the continental margin in the 
vicinity of the Barrow Current is a departure from the level bottom 
assumed. Theoretical studies of topographic effects on the baro¬ 
clinic instability problem cited in (IcVJilliams (1979) review show 
that the classical Eady waves studied here are not strongly influ¬ 
enced by depth variations and that instability may even be enhanced 
by such variations. It is reasonable that bottom topography would 
not seriously modify the waves studied here which are trapped pri¬ 
marily against the upper surface of the ocean. 
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ILLUSTRATIONS 


Fig. 1 - Profile of horizontal current through an Arctic Ocean 

eddy. Solid line is current speed and dashed line is direc¬ 
tion, 5 May 1975, 76°23'N 149°48'W (Manley, 1981)., 

Fig. 2 - Profile of temperature (T) , salinity (S) , and density (<3"t ) 
for the same eddy shown in Fig. 1. Solid lines are for condi¬ 
tions outside eddy in surrounding waters and dashed lines are 
for conditions inside eddy (Manley, 1981). 

Fig. 3 - Locations of arctic eddies observed during drift of four 
AIDJEX ice camps. Of the 98 eddies for which rotation sense 
was determinable, 95 were anticyclonic or clockwise. The 
three cyclonic eddies are circled. 

Fig. 4 - Phase speed and growth rate of baroclinic waves. H = 5. 
Solid line is for rigid upper and lower surfaces with no 
friction. Dashed line is for Ekman pumping at a free upper 
surface and rigid lower surface, W =0.05, Y =0.5. Dotted 
line is for Ekman pumping at upper and lower rigid surfaces, 

Y = 1. 

Fig. 5 - Phase speed and growth rate of baroclinic waves. H = 10. 
Symbols as in Fig. 4. 

Fig. 6 - Relative amplitude, M, and negative phase,— © , of fast¬ 
est-growing wave. Rigid upper and lower surfaces with no 
Ekman pumping. H = 5 and R = 0.8. 

Fig. 7 - Relative amplitude, M, and negative phase,—© , of fast¬ 
est-growing wave. Ekman pumping with free upper and rigid 
lower surface. 0=5 and £ = 0.7. 

Fig. 8 - Relative amplitude, M, and negative phase, — © , of fast¬ 
est-growing wave. Ekman pumping with rigid upper and lower 
surfaces. H = 5 and £ = 0.4. 

Fig. 9 - Geostrophic currents in the Arctic Ocean. Profiles (1) 

and (2) are in the central Arctic Ocean, (3) is near the shelf 
edge north of Alaska (Kusunoki, 1962). 

Fig. 10 - Profiles of VHisaia frequency in the central Arctic Ocean. 

Fig. 11 - Development of a vorticity profile with a subsurface max¬ 
imum from an initial exponential profile. Vorticity at the 
upper surface is kept at zero by frictional dissipation. 
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